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SUMMARY 



The equations relating the wing and tail loads are 
derived for a unit elevator displacement. These equations 
are then converted into a nondimens i onal form and general 
charts are given by which the wing- and tail-load-increment 
variation may "be determined under dynamic conditions for 
any type of elevator motion and for various degrees of 
airplane stability. In order to illustrate the use of 
the charts, several examples are included in which the 
wing and tail loadsare evaluated for a number of types of 
elevator motion. Methods are given for determining the 
necessary derivatives from results of wind-tunnel tests 
when such tests are available. 

INTRODUCTION 



Because airplane failures in which tail surfaces were 
involved have occurred in flight, recently, considerable 
impetus has been given to the tas.k of setting up more 
rational methods of evaluating tail loads. Particular 
interest has been shown in the analysis of dynamic tail 
loads associated with more or less sudden elevator motions. 

The problem of determining the dynamic tail loads in 
a. rational manner has been treated by many authors. Vari- 
ous approaches and assumptions have been employed, but the 
methods available at present are too lengthy to be suitable 
for the routine computations that would have to be made in 
design studies. This statement is particularly true if the 
critical types of elevator motion are to be varied consid- 
erably, from the simple types that have usually been 
treated. Although equations were given in reference 1 for 
determining the tail load with any variation of elevator 
motion, the equations were' not in the best form for making 
computations. It has been found recently, as a result of 
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a number of c ompixtat i ons , not only that the method of ref- 
erence 1 can "be shortened but also that some of the minor 
factors which were previously omitted can now be included 
in a method that will be suitable for use by designers. 

SYMBOLS . 

The following is a list of the symbols employed in 
this paper : 

V/ airplane: we ig.h.t , pounds 

g accelerat ion. of gravity, feet per second 3 
m airplane .mass , W/g, slugs 

S gross wing area including' area within fuselage,' 
square feet 

Sj. gross horizontal-tail area including that intercepted 
by fuselage, square feet 

b wing span, feet 

bj. tall span, feet . 

ky radius ..of gyr at i on- ab out pitching axis, feet 

■I". pitching moment of inert ia , slug-feet square 

length from center of gravity of airplane to aer'odyv- 
• namic center of tail* (negative for conventional 
.airplanes ) , f eet • 

V air plane velocity , f eet per second . 

p.. mass density of air, slugs per cubic foot 

q dynamic! pressure , pounds per square foot ^~pV' s ^ 

t tail, efficiency factor : (q t /q) 

L lift, pounds 

C L lift coefficient (L/qS) 

M moment, foot-pounds 
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C m pitching— moment coefficient of airplane without 
horizontal tail (Mb/qS^) 

a wing angle of attack, radians 

tail angle of attack, radians 

i t tail setting, radians 

6* elevator angle, radians 

€ flownwash angle, radians ( — a) 

y flight-path angle, radians 

9 angle of pitch (a + Y ) , radians 

K empirical constant denoting ratio of damping moment 
of c omple t e airplane to damping moment' of tail 
alone 

n airplane load factor 
t t ime , sec onds 

T aerodynamic time, unit = m/pSV 
p. airplane density ratio (-m/pSx^,) 

a,"b roots of tasic differential equation when they are 
imaginary 

m i , m 2 roots of basic differential equat ion when they 
are real 

K lf K 2 , Kz dimensional constants occurring in basic 
differential equation 

Kj 1 , K 2 ! , K3 1 nondimens i onal constants occurring in basic 
differential equation 

The notations a and a, 6 and b' , and so forth denote 
single and double differentiations with respect to. either 
t or T . 
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Subscr ipts 



0 


initial or selected value 


t 


tail 


max 


maximum value 


d 


down 




zero lift 


geo 


geometric 



THEORETICAL RELATIONS BETWEEN WING AND TAIL LOAD 



The mathematical treatment of the longitudinal motion 
of an airplane following an elevator displacement involves 
three simultaneous nonlinear differential equations . The 
correct analytical solution of these equations must, be ob- 
tained either by a series substitution or by step-by-step 
methods. A close appr oximat i sn to the correct solution. is 
obtained if it is assumed that, ili the. interval between 
the start of the maneuver and the attainment of maximum 
loads on the wing and tail surfaces, neither the initial 
velocity nor the initial attitude changes materially. 
These assumptions eliminate one of the three equations of 
motion and the trigonometric coefficients that occur in 
the other two equations. In addition, the assumptions 
agree with experimental flight results and have been gen- 
erally used in all treatments of the longitudinal motion 
of an airplane following a control deflection, 

If the sign conventions of figure 1 are used, the 
following equations will apply to the steady flight con- 
dition 

dC L 

W cos Y 0 - — ji a 0 qS = 0 (l) 

C m q— + i la 0 (l - Si j+ i t + --M 0 0"l t q')S t x t = 0 (2) 

b da^. L ^ ' do J 

Equation (l) represents the summation of the forces 
perpendicular to the instantaneous flight path and equa- 
tion (2) represents the moments about the center of grav- 
ity. 
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In accordance with the assumption that there is no 
loss in speed dtiring the pull-up, the corresponding dy- 
namic equations can be written as 

W cos(Y 0 + AY) (a 0 + Aa)qS--~i ri t qS t A6 + mYV = 0 (3) 

for the vertical forces. In this equation the term 
~™~ T lt ( l S t^ s is introduced to allow for the change in 

the Z force that will occur with elevator deflection. 
If the slope dC^/da is used for the complete airplane 

with the tail surfaces in place and elevator fixed, most, 
of the effect of the tail load on the vertical force will 
"be taken into account. 

The moment equation is 
V m da /*b da t I 0 V da J V da V Ji) 



+ i* + 



da t , .1 dC m t S t 2 



d8 



} ^ ill +• Of ♦ « 

(6 0 + A6 ) U t qS t x t : — -T) t q — A6 - mk y a e = 0 (4) 

J do D f 



In equation (4) the terra containing d is introduced to 
correct for the effect of time lag in downwash at the. 
tail, the term containing 6 is introduced to account 
for the change in tail angle due to rotation, and the 
dC mt S t a 

term ^ n^q — — A8 is introduced to account for the 

t 

moment due to elevator camber , Computations have indi- 
cated that in some cases it is necessary to include "both 
the camber term and the e le vat or -f or ce ' t erm . 

If equations (l) and (3) are subtracted from equa- 
tions (3) and (4), respectively, and if it is assumed 
that only a small change in attitude takes place (so that 
cos (Y 0 + AY ) S cos (Y ) ) , the following equations of mo- 
tion are obtained 

dC T d8 L 

mYV - — k Aaq.S — 1 n t qS t A6 = 0 (5) 

da do 
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*SS Aa, 5- ♦ i**Iia(l - !S\ - a f* 41 - | ft JL 
da b d a ; L \ da/ v da V vn^ 



da 
"dl 



t , 



&6 



dC m s t 2 ■ .. 

n t qS t X t " "d8 "t* -57 - mk y 2 6 = 0 (6) 



From figure 1 the following relations are seen to 
exist: 

8 = (a 0 + Aa) + (Y Q + AY ) 

B = a + i (7) 
6 = a + V 

Thus, from equations (5) and (?) 

♦ dC T S dCLf Sf 

. Y = 9 - a = -~— Aaq — + n t q £1 A6 (8) 

da mV d6 * mV 

and 

dC L . S dCL f S. . 

Y = 9 — a = = aq — + — -I r, q ~J 6 (9) 

da *mV d6 * mV 

If equations (8) and (9) are substituted into equa- 
tion (-6), the terms containing a, a, Aa,. A6 , and 6 are 
segregated and, if the resulting equation is divided by 
-inky 2 = -I, there is obtained 

L da| 2 I , P fc - da 2 mj * a [da q lb 

-da7 (T, t^ -r-'v 1 ~ ~ - i~ ^= 2 tA "l/irK^ ~i 

.'!!!i ( , t0 !t:_!!h iiVp q ] _ fe tq M (10) 

d8 * it t da t d6 2 -ml \| I dS * mV J 

The effect of the term containing 6 is small and may be 
omitted. Thus, Equation (10) can be written as 



a + K x a + K 2 Aa = K3 AS (11) 
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This equation is the equation for a damped oscillation 
with an impressed moment K3A6 where 



v pv r dC i.t s t x t 2 / 

Ll 2m Ua"" ~k~" H 



2*V ~ S 
da / da 



K 3 



2m Ida ky 3 b da t k Y s L\ da / da 2 m J J 

1 s t^t *t* ac h d % 5ni! p »t'»t' 1' . 

2m L d6 1 < k - 3 d& t b t k. y 2 da t d6 7^ 2 mk Y 2 J 



(Ha) 

The increment in wing l c fad, wing-load factor, and tail 
load can b$ found by solving equation (11) for ia and a 
by the usual methods, • The increment in wing load and wing— 
load factor could then be obtained from the equations 

dC L 

AL = J AaqS 

da 

do (12) 

An = t°-± 

da W/S 

It is seen from the bracketed term in equation (6) 
tha't,\in order to determine : the effective tail angle of 
attack, Aa t , at any time, the: pitching velocity and the 
rate of change of the wing angle of attack must first be 
known. If substitutions are made from equations (7) and 
(8) into this bracketed term, the increment in effective 
tail angl« of attack at any time is very closely given by 
the following equation 

L \ da da 2 m v /^ t J V \da Jr^J d8 • J 

The value of Aa^l given in equation (13) is to "be inserted 
in the aquation 

AL t = Aa t n t qS t (14) 

da t 

to ootain the tail-l?ad increment at any time. 
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Up to this point the equations and method are straight- 
forward and similar to the analysis previously presented in 
reference 1 with the exception that the rate of change of 
vertical force with elevator angle and the change ' in moment 
caused by tail camber have been- introduced into the equa- 
tions. These add i t i onal f ac t or s are usually small, hut they 
tend to gain in importance as the amount of static stabil- 
ity is increased. For the case of a very stable airplane 
their contributions may affect the results i.n the order of 
about 5 to 10 percent.' 

The solution of the differential equation of motion 
(equation (ll)) is not particularly difficult but would 
become rather tedious when the elevator motion is a com- 
plicated function of the time or when various types of 
elevator motion are to be considered. Also, in the form 
given, new computations would be required for each alti- 
tude and for each speed and the computations made for one 
airplane would not be applicable to another. 

The first difficulty can be avoided by evaluating the 
results for a unit instantaneous elevat or— angle change; 
then, since the equations are linear and the principle of 
superposition applies, Carson's or Duhamel f s integral 
theorem may be used (see reference 2 for application) to 
obtain results for any assumed elevator variation. The 
second difficulty can be partly overcome by selecting, as 
did Glauert, new units of time and length and presenting 
charts for the unit solutions of Aa and a for the 
various degrees of stability that would be obtained for 
center-of -gravity positions between the aerodynamic cen- 
ter and .the neutral point, 

In line with these ideas, the increment in elevator, 
angle will be taken as unity and the unit of time,- in- 
stead, of being taken as 1 second, will be taken as 

m 

T s seconds. The unit of length will be taken as 

pSV x+ • 

x t feet so that the unit of velocity will be or 

V/V where |a = — m/pSx t . Since x t is a negative quan- 
tity, with the system of axes used, p, will be a.posi- 
tive quantity the value of which may range from about 10 
to 100. 

Introducing the ab ove .quant it ie-s into equation (ll) 
allows a similar differential equation to be ~ob tained , 
which can be written as 
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a + Kj 'd + K a 'Ai = K 3 'A6 (1) (15) 
where A8(l) is a unit displacement and 

"dC L . S f Xt 3 / K de \ dC T "j 

K , 1 = Ij— ill _k n + __ )+ — L 

2t_da t S k y s 1 V"n*t da ' da -J 

2 Ida kv a . ' t s k a |\ da/ da Jn~ ~3 W 



* a k y 2 b da t S ky 2 LV da/ da ^ 

= ~ 2L "Is" n t ~ s ~ ~ "7~ r 't 

dC Lt AC Lt Knt s p x t 3 S t 2 



K3 . ..tti _ d !ii !t f tf . !S « t i_tl 

d6 fc S k y s d6 * b t Sk y 2 

Lt tf _P_ x t 3 S. 

da^ d8 Vtit" 2m Sk y : 



] 



(15a) 



The value of K l t is always positive and the value of 
Kg 1 is positive if the center of gravity lies ahead #f 
the rear neutral point. The rear neutral point is defined 
here as the position along the mean aerodynamic chord at 
which the center of gravity would have to be in order that 
the slope of the moment curve for the aaittftlfiAe. .airplane 
at out this point be 0. The value of dC m /da that is 
used in this report is taken a "bout a forward neutral point 
(with tail off), which has teen called the aerodynamic 
center. The quantity Kjj 1 is always negative and depends 
only on the geometric and aerodynamic qualities of the 
tail. 



The solution of equation (15) can" take any one of 
three forms, depending on whether both roots of the auxil- 
iary equation are real and unequal (ra lt m 2 ), real and 
equal (n^ = to 2 ), or imaginary in the form of a+ih. With 
the stipulation .that the canter of gravity be forward of 
the rear neutral point, the motion indicated "by equation 
(15) always subsides and the solutions for Aa and d 
are as f ollows : 
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Unequal real roots m l £ 



m. 



a * K 3 ig (1) 



e 



K 3 f 




I \2 



— — sinh 



^-__V K - It 



2 y 



k 3 ' - 



+ ^("s" ) " Ks ' cosh 



K 3 '8(l) 



s inh 



'( ___ > - IC - ' 




^ (16) 



aqual real roots m l = m 2 
K 3 '6 (1) 



1 - e 



K 3 "6 (1) 



a = 



K 3 ' 



l ' N , t e 



+ i 



(16a) 



Imaginary roots a±ib 



K 3 '5 (1) 
Aa = < 



-.lilt 



1 + e 



cos I / K s ' 



3 7K,.-(5|-j 

_____ 



\ 2 



sin ./ _ fl « 



a = ______2__. 



sin 



-(¥) 



(16b) 
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In these solutions, the boundary conditions are at 
T = 0, Aa = a = 0, A5(l) = 1.0. References to equations 
(5) and (7), however, indicate that the "boundary condi- 
tions should he at T = 0, Aa = 0, A8(l) = 1.0, e = 0, so 
dC T S t T) t 

that a = -Y = — - ~ — . The inclusion of these condi- 

d6 S 3 

tions complicates the solution and introduces factors that 
prevent the presentation of results in a few "basic charts. 
Actual plots of the unit-solution curves obtained with 
either boundary condition indicate, in examples that have 
been tried, such small differences that the two curves 
can he hardly distinguished. For these' reasons , the sim- 
ple boundary conditions have been used. 

It has been found by direct substitution that the. 
value of 'K 1 % will range from, about- 5 to 9 in the case 
of conventional airplanes. (See equation (15a).) Similar 
substitutions for- K 2 1 indicate that this quantity may 
range from about 2 to about 300 when all possible values, 
of p. and p are considered. There are, however, com- 
pensating factors that enter into the problem so that the 
likely range of K 2 r is much smaller than this even when 
the possible present-day extremes of the separate items 
are c ons ider ed , 



CHARTS F OR DETERMINING Aa AND a 



Charts are given in figures 2 to 6 showing the varia- 
tion of AaXs'/Ks 1 and a/K 3 f against aerodynamic time 
T for all values of K 1 1 and K 2 ! that are likely to 
occur. The charts given apply as long as K 2 1 remains a 
positive quantity, which will always be the c?se when 
there is a small margin of static, stab i lity , namely, when 
the center of gravity is ahead of the rear neutral point. 
According to the bracketed term of equation (15a), the 
center of gravity could be slightly behind the neutral 
point and the motion given by equation (15) would still 
subside because of the . greater stability which the air- 
plane has on a curved path. 



USE OF THE CHARTS IN A TYPICAL EXAMPLE 



In , order to illustrate the generality of the charts 
given in figures 2 to 6, an example is worked for a 
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typical fighter airplane which is now under investigation 
for tail loads. The necessary geometric and aerodynamic 
characteristics of this airplane are as follows: 



Geometric 

Gross wing area / S t sq ft . 300 

Gross horizontal tail area, S^, sq ft. 60 

Airplane weight, W , lb ..... . 12,000 

Wing span, 1), ft 41 

Tail span , b t , ft 16 

Radius of gyration, ky , ft " 6.4 

Distance from aerodynamic center of airplane less 

tail to aerodynamic center of tail, x t , ft . . .-21.0 

Aerodynamic 

,dC ; 

Slope of airplane lift curve, — - , radians 4.87 

ac Lt da 

Slope. of tail lift curve, i, radians 3.15 

da t 

Down wash f act or , d e/ da . , . , . . . . 0.54 

Tail efficiency factor (q^/q),Ti t . . . . 1.00 

Empirical airplane damping factor, K ...... . . 1.1 

Elevator effectiveness factor, - f radians .... 1.89 

d8 

Rate of change of • t ai 1 moment with camber due t o. 

AC m r 

elevator angle, -g^ — , radian . . . 0.57 

t 

Rate of change of moment coefficient with angle of 
at tack . 

(a) center of gravity, 30 percent, radian . . . 0.703 

(b) cen.ter of gravity, 25 percent, radian , . . 0.475 

It was determined from tests of this airplane that 
for the conditions desired the slope of the moment curve 
per radian for the. airplane less tail could be given by 

= _o . 665 + 0. 0445 c . g.. 

da 

Substitution of the geometric and aerodynamic values 
into equa-tion (15a) and the assumption that results are 
required for an altitude of 19,100 feet (p = 0.0013.06) 
give the following values for K s f , and K 3 1 : 



13 



K 1 ' = 8.0 

K 2 l (e.g. at 30 percent) = 20.0 
K s ; -(c.g. at 25 percent) = 40,0 
K 3 ' = -100. 0 

For these values of K x 1 and K 2 1 , the variation of 
K ^ 1 d 

Aa — — and with T can "be obtained from figures 

K 3 ! K 3 ' 

5(a) and 5(b) for an instantaneous unit elevator deflec- 
tion... A slight amount of labor can be saved at this 
stage if the curves are taken directly from these figures 
onto a work sheet (see middle group of curves in fig. 7) 
without transforming them into curves of Aa and d. 
The transformation can, of course,- be accomplished im- 
mediately by multiplying the ordinates of the curves ob- 

tained by and K 3 1 , respectively. It has been 

&2 1 

found more convenient, however, to make the change— over 
as a final step. 

The next step in the procedi—e is to plot the assumed 
elevator-motion curve on the work sheet using the same 
abscissa (t). This change is accomplished by dividing the 
actual assumed time variation of elevator deflection by 
the factor m/pSV in order to obtain the variation in 
aerodynamic units. For an indicated speed of 400 miles 
per hour at 19,100 feet, the factor m/pSV for the 
airplane in question would be 



.V//S • 40 



(gp)V</p 0 /p (88/60) 0 . 0420X400X1.349X1 .466 



= 1.202 seconds 



The determination of' Aa — ^— and — — at any time 

T due to the assumed elevator motion is then found by 
the following graphical construction. This construction 
is essentially that given in reference 2 except for minor 
modifications that were found to be worth while in effect— 

K * a 

ing the computations. The values of Aa — - and at 

the aerodynamic time t Q = 1, for example, due to the 
assumed elevator motion is found' as follows: 
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1. First, the point on the A6 curve at T = 1 is 
projected horizontally in "both directions until it strikes 
the 45° lines. The intersections with these lines are 
then projected or deflected vertically until they inter- 
sect the -hor i zont al projections of the values of d/Ks 1 

and Aa — — at the time T equal to zero. The- points 

labeled ~(T) are thus established. 

2. The ordinate of the A6 curve at, say, T = 0.8 
is next projected horizontally as before until it strikes 
the 45°- lines, where it is reflected and points of inter- 
section- with the horizontal projections of the values . of 

Aa . and a/K^ 1 at T = 0.2 . second are established. 

^ K3 1 



The point s labeled (2) are obtained in this manner. 
Other points labeled (5) 1 0 (5) for t = 0.2, 0.1, and 0 
are then similarly obtained to complete the curve for the 
example chosen and for the time t 0 = 1.0. Note that the 
addition of ■. T on the elevator curve and t on the unit- 
function curve always equals T 0 , Curves are then drawn 
through these' points and the areas under them are propor- 
tional. to. Aa and a for the aerodynamic time of 1 unit. 

: ; .3 , /The -areas ar e found by • integrating .in the direc- 
tion shown, .that is , .© , ©,©.,. ® , .and so forth. 
It ie .. important .to -follow .in this direction in order that 
negative, areas, if ' they .should . occur , may 'be properly 
taken into account.. If .. the " point s are followed in the 
order noted and a counterclockwise path is followed' in 
enclos ing ' the area, the val ae is pos itive ro^ai'dln gr, of 
the quadrants involved and vice versa for clockwise in- 
tegration. . When a figure-of-eight area is involved, the 
same statement also applies. The areas are then converted 
to Aa and & by multiplying the number of square units 
(square inches or square centimeters) by the appropriate 
conversion factor's, which are . obtained . by multiplying the 
ordinate scales of Aa or - a by the ordinate scales of 
AS curves as the case may be. 

. 1 4. Other curves Are similarly drawn . in for the dif- 
ferent time intervals, T Q , at which the values of Aa, 
and d are desired. For example, see the heavy curve 
drawn in for the time t = 2.2 with points labeled CD, 



ijLl , and so f or th . 
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5. After a sufficient number of time intervals are /" 
considered and the resulting areas determined, the final 
step is to convert the areas into the values of Aa and 
d associated with the elevator motion assumed. The var- 
iation of load factor and tail— load increment may then 
"be found by substituting for Aa and & in.to equations 
(12), (13), and (14). It is convenient at this stage to 
arrange the results in tabular form and to convert from 
timeTtotimet. 

Figure 8 gives the final variation of the load fac- 
tor and tail-load increments for the example of figure 7. 
The type of elevator motion used in these figures is that 
previously suggested in reference 1 in which the horizon- 
tal tail would be designed. to withstand a maneuver in 
which the V— G diagram. from maximum' posit ive to maximum 
negative g would be covered. The duration of the time 
intervals at which the elevator was- held at maximum val- 
ues iASmax was adjusted so that the full acceleration 
corresponding to each .elevator -throw would be reached. 
The rates of movement were purposely t?ken quite high in 
order to obtain as large downtail or control loads as 
possible without exceeding the speed at which the pilot 
Might move the controls. The relation between the ele- 
vator throw A6 filax and the load factor increment An 
that is finally reached is given by 

K 2 ' W/S 

A6 inax = An 

— 

' da 

This relation is easily obtained by substituting the val- 
ues for , Aa from equations (L"6) (with t large) into equa- 
tion^). 

In order to obtain this range in acceleration, it is 
not. necessary that the pilot restrict himself to the type 
of elevator motion assumed, as lie may actually move the 
elevator twice as far as is necessary and check the motion 
earlier so as not to overshoot the desired acceleration, 
Such a motion is illustrated by the results given in fig- 
ures 9 and 10 for the same airplane (K^ = 8 , K 2 1 = 20) 
but' with the elevator motion required to cover approxi- 
mately the same acceleration range. 

It will be -noted in figure 10 that the maximum accel- 
eration reached with the elevator motion assumed was 8.75g 
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instead of -8g. A number • of tr ial computations made in 
connection- with the " f igure • indicated- that , other things 
being equal, a delay of as little as 0,06 second in the 
time of- the elevator reversal would cause the accelera- 
tion to. overshoot "by 1.5g.' This' delay indicates that 
the par t icular . type of elevator motion shown in figure 
10 would probably never be used by a pilot in a high g 
pu-ll-out and where the elevator motion is small because 
of the extremely fine timing requiredto prevent over- 
loading . 

Comparison of the results of- figure 10 with those 
of figure 3 for the same values of K s 1 , indicates a much 
more rapid variation in load factor for the type of mo- 
tion used in figure 10, Jones and Fehlner , in reference 
3, have shown that the transient effects of wing wake 
on the tail are likely to be severe only when the rate of 
change of wing angle of attack is great . These effects 
are not included in the method given because in the usual 
case they- apparently are of little importance in the de- 
termination of the cr it ical— maneuver tail loadi In order 
to illustrate the combined effects of aerodynamic lag and 
transient wing wake on tail loads, a portion of the tail- 
load curve in the critical region, including transient 
effects, .has been computed by R. T. Jones for the case 
illustrated in figure 10. The comparison is given in fig- 
ure 11 where it will be seen that even in this particu- 
larly severe case the discrepancy amounts to only about 
10 percent on the important maximum loads. 

Because actual elevator motions are almost certain 
to be less severe than the one illustrated, the transient . 
effect will be less than that shown and within the limits 
of accuracy with which some items entering into the. .com- 
putations are known. Per this reason, and because 'of"«the 
increased mathematical complexity that is introduced by 
its inclusion, transient effects are omitted. This state- 
ment ,. however , cannot be assumed to apply to the gust,-., 
condition wherein the angle-of-at tack changes may occur, 
jnore' rapidly. ' m 

POSSIBLE SHORTCUTS ... 

The preceding section illustrated a general proce- 
dure that can be' f ollowed where one or. two el-evator motions 
are to be inve s t igat ed^but ; there are a number of variations 
which might have' been used. If, however, the effects of a 
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fairly large number of elevator motions are to "be inves- 
tigated at a given speed and altitude, the following 
method can be used with a saving in time. 

1, Determine a unit tail load AL^ . and wing load 
An directly by using values of" a and Aa obtained 
from figures 2 to 6 and substituting these values into 
equat i ons (12 ) and (14 ) . 

2, Plot the values An and AL^. of step 1 as a 
function of time instead of plotting Aa and d as a 
function of time as in figures 7 and 9. 

3, Plot the elevator-motion curve to be investigated 
and proceed as before. The areas now obtained will give 
directly the increments ' of tail load and wing load. 

A quicker me t hod of obtaining the areas , in s ome 
cases, and one that is readily apparent after a little 
experience is gained, is to prepare grids for e valuat ing 
either Aa and & or An and AL^ . The abscissas of 
the vertical lines of the grids are simply the ordinates 
of the "considered e le vat or— mot ion curve taken every 0.1 
or 0.2 second, say; the ordinates of the horizontal lines 
of t he" gr id are then the ordinates of the respective unit 
curves. Figure 12 shows such a : f ;;rid for determining 

Aa — — for the conditions of figure 7(a). In order to 

obtain the value of Aa — — for the specific time of 

t 0 = 1.0, the points of intersection of the horizontal 
and vertical gr id lines adding up to 1.0 ar e c onne c t ed . • 
Such a curve is shown in figure 12 for comparison with 

the similar one given in figure 7(a). 

1 

DETERMINATION OF THE NECESSARY AERODYNAMIC DERIVATIVES 

The accuracy, with which the load increments may be 
determined for a given elevator motion depends largely 
upon. .the accuracy with which .certain aerodynamic charac- 
teristics are known and, in some measure, on how well 
these char act er is t i c s may be appr oximat ed by a s tr ai ght . 
line. The values required for the computation may be 
obtained with sufficient accuracy from wind-tunnel tests 
of a model in which the lift, drag,' and moment are meas- 
ured with and without the tail in place and with the 
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power condition for which calculations are" to "be made. 
Lift and moment measurements should also "be made for a 
range of elevator angles from ±10°; then, with elevator 
fixed, the lift and moment variation v/ith tail setting 
should be determined through a range of about 5°. 

The value of dCx,/da to.be used should be that . 
obtained with the tail in place and should be based on 
the gross wing area. By the use of this value, most of 
the effect of tail load on normal acceleration will be 
taken into account. 

The value of dC m /da to be used is the slope of the 
moment curve with the. horizontal t a il . r'e move d . Usually, 
the moment variation is taken with respect to a selected 
cent er~of ~-gr avi t y position, but it is desirable that the 
variation of dC m /da with c ent er-of — gr av ity position be 
determined for at least two center~of— gravity positions. 
For the conventional fighter airplane it appears that the 
critical total down load at the tail will occur with the 
center of gravity in its most forward position during 
dive pul l~out s at high al t i t ude and at . t he 1 imi tingle peed-. 
The maximum up— load at the tail is likely to occur during 
pull-ups from high-speed level flight with the center of 
gravity in its most rearward position and at ,a relatively ' 
1 ow alt i tude , . 

dCL t 

The value .of the factor r\, can be obtained 

t <*a t 

from moment differences obtained at the same angle, of 
attack from two settings i^ of the horizontal tail 
plane . Thus , 

' tda t Ai t bS t x t 

It will not generally be necessary to separate the factor 
71^ but this separation could be accomplished by reference 
to tests- of isolated tail surfaces of a similar plan form, 
Reference 4 gives results for tests of a number of isolated 
tail surfaces . 

dC L 

The previous va lue ob ta ine d for ti , — — — can be used 

dC Lf 

to find the elevator effectiveness da^/dS or from 

•either the moment or the lift differences obtained from 
tests in which the elevator angle was varied. 
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Moment difference's 



flat = ^Cm av S ; 



d6 dC Lf 

n. S t x t bA6 

* da t t * 



dC Lf AC sS 



11 1 ~d6 



S t x t bA8 



Lift differences 



1*1 _ _1 '^^LajrS 



d8 dc Lt 



It 



— S+A6 

da t * 



AC T S 



v 



T). £ = — 

t d8 S t A5 



The differences AC L and AC m are to be taken at the 
same angle of attack. In general, the moment differences 
will prove to "be the most reliable because the quantities 
involved are larger. 

Similarly, the dcwnwash factor can be determined on 
the basis of either moment or lift differences with and 
without the tail in place together with the previous value 

dC L 

obtained for n, £ 

t da t 

Moment differences 



/\ deA ». /dC m > S 2 1 

[ 1 - — j = A • — - \ ± 

V da / \da / S t "bx t dC L 

r, . i 



t 



da t 



where A( — - ) is the increment in the slop 

\ da / / da 

caused by the addition of the tail. 
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Lift differences 



d£ \ S [" f6.C k \ 

da/ " dC L t L V da J 

S t n t Ta~" 



tail on 




tail off J 



The values obtained from the moment differences are the 
more reliable . 
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Figure //*- Effect of transient wing woke on The tail load. 
(Measure rv/'rh <4Ctl' scale) 



